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Abstract 

Systems with bond disorder are defined through lattice Hamiltonians that are 
of pure nearest neighbour hopping type, i.e. do not contain on-site contribu- 
tions. They stand representative for the general family of disordered systems 
with chiral symmetries. Application of the Dorokhov-Mello-Pereyra-Kumar 
transfer matrix technique [P. W. Brouwer et al, Phys. Rev. Lett 81, 862 
(1998); Phys. Rev. Lett. 84, 1913 (2000)] has shown that both spectral 
and transport properties of quasi one-dimensional systems belonging to this 
category are highly unusual. Most notably, regimes with absence of exponen- 
tial Anderson localization are observed, the single particle density of states 
exhibits singular structure in the vicinity of the band centre, and the man- 
ifestation of these phenomena depends in an apparently topological manner 
on the even- or oddness of the channel number. In this paper we re-consider 
the problem from the complementary perspective of the non-linear fi-model. 
Relying on the standard analogy between one-dimensional statistical field the- 
ories and zero-dimensional quantum mechanics, we will relate the problem to 
the behaviour of a quantum point particle subject to an Aharonov-Bohm flux. 
We will build on this analogy to re-derive earlier DMPK results, identify a new 
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class of even/odd staggering phenomena (now dependent on the total number 
of sites in the system) and trace back the anomalous behaviour of the bond 
disordered system to a simple physical mechanism, viz. the flux periodicity of 
the quantum Aharonov-Bohm system. We will also touch upon connections 
to the low energy physics of other lattice systems, notably disordered chiral 
systems in and 2 dimensions and antiferromagnetic spin chains. 

PACS numbers: 72.15.Rn, 72.10.Bg, ll.30.Rd 
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I. INTRODUCTION 



In 1979, Wegner and Oppermann introduced a variant of the disordered lattice An- 
derson model they termed 'sublattice system'. The sublattice system differs from the generic 
Anderson model in that its Hamiltonian does not contain on-site matrix elements, i.e. is 
of pure hopping type. For a long time this species of disordered electronic systems went 
largely unnoticed. The status rapidly changed PHl3| after two aspects became generally 



appreciated: first, models with sublattice structure occur in a number of physical appli 
cations. The random flux model, lattice QCD models [|14|, random antiferromagnets [jl5[ 



models of gapless semiconductors [T^ and effective models of transport in manganese oxides 
ir?! are of sublattice type or at least acquire sublattice structure in limiting cases. Second, 
and contrary to naive expectations based on universality and insensitivity to details of the 
microscopic implementation of disorder, the low energy properties of the sublattice system 
differ drastically from those of the generic Anderson model: 

• In contrast to the Anderson model, average spectral and transport properties of the 
sublattice system sensitively depend on the value of the Fermi energy, Ep. For Fermi 
energies far away from the centre of the tight binding band, Ep = 0, the sublattice 
system falls into the universality class of standard disordered electron systems (which 
simply follows from the fact that a Fermi energy Ep 7^ can be interpreted as a 
constant non- vanishing on-site contribution to the Hamiltonian.) However, in the 
vicinity of zero energy drastic deviations from standard behaviour occur: 

• In dimensions d < 2, the average density of states (DoS), i^{E), exhibits singular 
behaviour upon E approaching the band centre. 

• Perturbative one-loop RG calculations as well as the analysis of Ref. indicate 
that right in the middle of the band the 2d system is metallic, i.e. does not drive 
towards an exponentially localized phase. 

• Several phenomena of apparent topological origin are observed. E.g., the DoS proflle 



of sublattice quantum dots (systems in the 'zero-dimensional' limit) sensitively de- 
pends on the total number of sites of the host lattice being even or odd. Similarly, 
the properties of quasi one-dimensional sublattice systems depend on the number of 
channels, A^, being even or odd. For even the transport behaviour is conventional 

- conductance exponentially decreasing on the scale of a certain localization length ^ 

- whereas the energy dependent DoS vanishes in a close to linear fashion. In contrast, 
for N odd, the wire at = is metallic, i.e. the length-dependent conductance decays 
only algebraically. At the same time, the DoS diverges upon approaching the band 
centre. 

A schematic summary of the band centre phenomenology of sublattice models is displayed 
in table |I|. 

All these phenomena root in the fact that the sublattice Hamiltonian possesses a dis- 
crete symmetry not present in the Anderson model: [a3,H]+ = 0, where [ , ]+ is the 
anti-commutator and a site-diagonal operator that takes values +1/ — 1 on alternating 
sites. The presence of this 'chiral' symmetry implies that sublattice systems fall into sym- 
metry classes different from the three standard Wigner-Dyson classes 'unitary', 'orthogonal' 
or 'symplectic'. To be specific, let us focus on the simplest case of a sublattice system 
with broken time reversal invariance but good spin rotation symmetry (the analogue of the 
Anderson model of unitary symmetry.) Chiral systems fulfilling these two extra symme- 



try criteria belong to a symmetry class denoted Alll in the terminology of Ref. [18|. An 



alternative denotation, coined in the paper [2C], is ChGUE for 'Chiral CUE'. 
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FIG. 1. Schematic overview over phenomenology of sublattice systems. N : number of channels 
of a quasi one-dimensional quantum wire, L : number of sites of a sublattice quantum dot. 



As with conventional disordered electronic systems of Wigner-Dyson symmetry, universal 
transport and thermodynamic properties of systems of class Alll can be described in terms 
of effective low energy theories. E.g., the results for sublattice quantum wires summarized 
above have been obtained within a symmetry adapted variant of the Dorokhov-Mello-Kumar- 
Pereyra (DMPK) transfer matrix approach [^|TT],|TU| . This theory differs from the standard 
cases of unitary symmetry in that the transfer matrices describing the propagation through 
the system take values in a different target space. Similarly, the general field theory approach 
to disordered electronic systems, the nonlinear cr-model, has been extended to systems of 
class Alll [3,^p!3l,|6| , too. Like its conventional relatives, the Alll variant of the cr-model 
is a matrix field theory whose internal structure depends on the specific implementation 
(boson replicas, fermion replicas or supersymmetry.) Previous studies of these models have 
focused on the two-dimensional case MM, where information on long range behaviour can 



be obtained from renormalization group calculations, or on the zero-dimensional case |21 



where the model can be evaluated rigorously by full integration over the zero-mode manifold. 

It is the purpose of the present paper to analyse the intermediate, quasi one-dimensional 
variant of the field theory. As mentioned above key aspects of the phenomenology of quasi 
one-dimensional sublattice quantum wires have been discussed previously within the frame- 
work of the DMPK approach, and it is near at hand to ask what motivates revisiting the 



problem. Referring for a more substantial discussion to section [I A| below, we here merely 
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mention a lose collection of points. The field theory approach enables one to approach the 
problem from a comparatively broad perspective. Specifically, the one-dimensional variant 
of the model is but a representative of a larger family of 'chiral' a-models. This makes pos- 
sible to relate the behaviour of the one-dimensional system to the extensively studied zero 
and two dimensional cases. Further, the Green function oriented cr-model formalism enables 
one to 'microscopically' implement coupling operators connecting the wire to external leads. 
(Within previous DMPK formulations, the coupling has been treated somewhat implicitly, 
see however |^2|.) Unexpectedly, we will observe strong, albeit universal sensitivity to the 
modelling of the coupling and yet another class of staggering phenomena. The origin of 
these effects, and their connection to the channel number staggering mentioned before will 
be discussed below. Finally, the cx-model of the time reversal non-invariant sublattice is the 
by far most simple of all ten [|1^] nonlinear cr-models of disordered systems: it has only four 
degrees of freedom, two Grassmann and two ordinary integration variables, the minimal set 
needed to construct a supersymmetric matrix model. This makes it an ideal tutorial system 
on which generic properties of the field theory approach to disordered quantum wires can 
be studied. We have tried to pedagogically expose several of these aspects, particularly the 
analogy one-dimensional field theories vs. point-particle quantum mechanics which plays an 
all important role in the present context. Nonetheless, the analysis below will be at times 
technical. To make its results and various qualitative connections generally accessible, the 
following section contains a synopsis of the paper. 



A. Summary of Results and Qualitative Discussion 

Consider the system depicted in figure Fig. |[ a sequence of L sites (alternatingly desig- 
nated by crosses and circles) each of which supports electronic states, or orbitals (repre- 
sented by the vertical stacks of ovals.) Nearest neighbour hopping is controlled by a regular 
tight binding contribution, diagonal in the orbital index, (the horizontal line segments) plus 
some bond randomness that connects different orbitals (the hatched areas). As in Refs. [Q, 
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we allow for some 'staggering' in the tight binding amplitudes, i.e. the hopping amplitudes 
regularly alternate in strength (the alternating bond length.) At either end, a number of 
sites is coupled through some tunnelling barriers (horizontal hatched areas) to leads. 



2 



2 
1 



w 

m- 



^ — O; (X 



L-l 



L 



4^ 



o) :x — o; (X — et 



\/\/\/\/ 



\/\/\/\/ 
\/\/\/v 



FIG. 2. Quasi one dimensional sublattice system. The system is realized in terms of a sequence 
of sites, alternatingly designated by crosses and circles, each of which carries N states (the ovals). 
Hopping between the sites is mediated by a regular Hamiltonian (the lines connecting the ovals) 
plus a weak random contribution (hatched areas). At either end, a number of sites is coupled to 
continuum states (the wavy lines) propagating in leads. 



In this paper, three different regimes will be considered: (i) the 'quantum dot' case, 
defined through the criterion that the time to diffusively propagate through the system is 
shorter than the Golden rule escape time into the leads, (ii) a diffusive regime, where the 
order of these time scales is inverted but localization effects do not yet play a role, and (iii) 
the regime of long systems with pronounced Anderson localization. 

Beginning with the quantum dot case (i), we find that both transport and spectral 
properties sensitively depend on the number of sites L being even or odd. Specifically, the 
DoS exhibits a singular spike in the band centre if the number of sites is odd (and the coupling 
to the leads switched off). Away from zero energy, v{e) is strongly suppressed up to values 
u NA, where A is the mean level spacing. The conductance equals 'yN/2, where 7 is a 
measure for the strength of the lead coupling, as for conventional quantum dots. In contrast, 
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for an even number of sites, the DoS shows the spectral 'microgap' of width A characteristic 
for finite random systems with chiral symmetries [|T^. Curiously, the conductance vanishes 
provided that only one site at either end is connected to the leads. The strong sensitivity 
to the total number of sites in the system disappears if more than one site in the terminal 
regions is coupled to leads (cf. Fig. 0.) Summarizing, the tendency of sublattice systems 
to exhibit staggering phenomena pertains to the zero- dimensional limit. However, unlike in 
the localized regime, the relevant integer control parameter is the number of sites L and not 
the number of channels A^. A qualitative interpretation of these phenomena will be given 
below. 

In this paper, only limited attention will be payed to the intermediate diffusive regime 
(ii). It is likely that the staggering phenomena observed in the zero-dimensional regime will 
have interesting, albeit non-universal extensions into the diffusive regime. We here avoid 
the confrontation with these effects by connecting several sites to the leads. This leads to 
equilibrated behaviour with Ohmic conductance, as for ordinary wires. As for the density 
of states, we have explored the influence of spatially fluctuating diffusion modes on the 
spectral microgaps discussed above. (Semiclassically, the spectral gaps observed in chiral or 
superconductor systems can be interpreted as due to an accumulation of diffusion modes (see 
Ref. 1^ for a detailed discussion of this point). For large energies e ^ A, these modes can be 
treated perturbatively by standard diagrammatic methods.) Surprisingly, no corrections are 
found up to and including three loop order. This is a speciality of the one-dimensional case. 
For a two dimensional system diffusive modes would lead to a modulation of the spectrum 
on the scale of the Thouless energy. 

In the localized regime we reproduce the results found earlier within the DMPK approach: 
for an even number of channels, the conductance decays exponentially on the scale of a 
localization length o(. Nl, where / is the mean free path. In contrast, for an odd number 
of channels algebraic scaling, g ~ L^^/^, supported by one delocalized mode is observed 
A comparably strong even/odd dependence is observed in the behaviour of the DoS. 
For A^ even, the DoS vanishes at zero energy as z/(e) ~ — |s| ln(|s|), where s = vre/A^, and 
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the characteristic scale of the gap, A^, is the level spacing of a system of length L ~ ^. 
This behaviour has a simple interpretation: roughly speaking, a system with L ^ ^ can 
be viewed as a sequence of L/C, decoupled 'localization volumes'. On small energy scales 
e < A^, dynamics within each of these is approximately ergodic. One would then expect 
the DoS to be gapped, as for the sublattice quantum dot, with A^ as the characteristic level 
spacing. Generalizing an earlier result ||23|j2^ for the specific case = 1, it has been found 



that for an odd number of channels the DoS diverges as z/(e) ~ l/{\s\W\s\) ini. This 



accumulation of spectral weight can be interpreted in two different directions. First, it is 
near at hand to view the algebraically decaying conductance observed in the odd case as 
a resonant tunnelling phenomenon: the principal tendency to localize is outweighed by the 
high density of states in the vicinity of zero energy. (This picture was first suggested by 
V. Kravtsov.) Second, it is tempting to interpret the zero energy peak as an generalization 
of the singular spike found in the L-odd zero-dimensional case. Unfortunately, we are not 
aware of a qualitative picture explaining this analogy. At least technically, however, both 
phenomena can be traced back to a common origin. Finally, a periodic modulation, or 
staggering, of the hopping amplitudes can be employed to continuously interpolate between 
the iV-even and A^-odd case, respectively. 

We next briefly outline the field theory route to exploring the above phenomenology. 
Within the fieldtheoretical approach, long range properties of the system are described by 
a functional integral with action 



rL 

S = I dr 







str {OTdT-^) - str {TdT-^) - i- str (T + T'^] 



5'Gade + S^. (1) 



Here we have introduced a continuum variable r G [0, L] replacing the formerly discrete site 
index, / is a parameter related to the staggering strength, s = €7tu, where u is the bulk 
metallic density of states. Further, T is a matrix field taking values in GL(1|1), i.e. the group 
of two-dimensional invertible supermatrices, and 'str' is the supertrace. Finally, ScadefT"] = 
C str ^(T9T^^), where C is some small constant and St is a contribution describing the 
coupling to the leads. 



The action S[T] defines a one-dimensional representative of the general family of 'chiral' 
non-linear cr-models. In contrast to its well investigated zero- and two-dimensional relatives 
- much of the results summarized in table ^ have been derived within these models - the Id 
variant has not been explored so far. 

Much of our analysis of this model will rely on the standard analogy between one- 
dimensional statistical field theory and zero- dimensional quantum mechanics: S[T] can be 
interpreted as the quantummechanical action of a point particle propagating on the super- 
symmetric target manifold, in our case GL(1|1). The first term of the action represents a 
kinetic energy, the third term a potential, and the second term, linear in the 'velocity', cou- 
pling to a constant vector potential of strength (A^ + /)/2. Quantum analogies of this type 
have been proven a powerful technical tool in previous analyses of the standard cr-models 



| 25| . However, the present case is special in that the target manifold is so simple that an 

intuitive interpretation of the quantum picture becomes straightforward. Indeed, the fields 

(u p\ 

T have the explicit matrix representation T = , where p, a{u, v) are Grassmann 

\a v) 

(commuting) variables, to be compared with the four or eight dimensional matrices entering 
the standard a-models. Later on we will see that convergence criteria constrain the compo- 
nent u to be positive real, while v = exp(i0) is a pure phase. Thus, temporarily leaving the 
Grassmann variables aside, our model describes quantum propagation on a (half)line and 
on a circle, respectively. Notice that the latter component is topologically non-trivial. 

At this stage, the role of the vector potential contribution to the action becomes evi- 
dent. While inessential in the non-compact sector of the theory, in the compact, circular 
sector, it desribes the presence of an Aharonov-Bohm type magnetic fiux. This analogy 
explains the presence of phenomenona periodic in the channel number. An even number 
of channels translates to an integer number of fiux quanta through the ring, which has no 
effect. However for A^ odd or, alternatively, a finite staggering parameter /, a fractional 
fiux pierces the system and this influences both, spectrum and dynamics of the quantum 
system. To develop the picture somewhat further, notice that the conductance, essentially 
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the transition probability through the system, maps onto the Green function of the quantum 
system evaluated at imaginary 'time' L. Imagining the latter represented through a spectral 
decomposition, the large L behaviour depends on the low lying portions of the quantum 
spectrum, in particular the discrete, and flux periodic level structure of the compact sector 
of the theory. Later on we will see that for half integer flux (i.e. odd) there is a zero-energy 
level (— s> absence of localization) while for all other magnetic configurations the spectrum is 
gapped exponential localization.) This mechanism, and its relevance for the localization 
behaviour of the system was first analysed by Martin Zirnbauer |2^. Similar but slightly 



more elaborate arguments can be used to understand the profile of the DoS. 

We finally mention some intriguing parallels to the physics of the antiferromagnetic spin 
chain. According to Haldanes conjecture, a chain of spins with half-integer (integer) S is 
in a long range ordered (disordered) phase It has also been found (see Ref. for 



review) that for a chain with staggered hopping amplitude j, the strict integer /half-integer 
pattern is violated, e.g. an integer chain can be fine-tuned to an ordered phase. Technically, 
the system is described by a cr-model with a topological term not dissimilar to the one 
above. The differences are that in the spin case (a) the base manifold is 1 + 1-dimensional 
(a quantum chain), (b) the target manifold is the two-sphere ~ SU(2)/U(1), and (c) the 
topological term classifies field configurations according to the number of coverings of the 
sphere (instead of winding numbers around the circle, as in our case). In the spin case, 
the coupling constant of the topological term is given by S + j, i.e. spin replaces channel 
number and staggering plays a similar role as in our case. (In fact, the linear coupling of 
the topological term to the staggering amplitude follows, independently of the model, from 
parity-conservation arguments to be discussed below.) Beyond these apparent technical 
parallels, the connection between the quenched disordered sublattice and the spin chain, 
respectively, is not understood. 

The rest of the paper is organized as follows. In section |I| we quantify the definition 
of the model. Its field theory representation is introduced in section |T|. In section ^ the 
cr-model transfer matrix approach, i.e. the representation of observables in terms of the 
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quantum Green function is discussed. This formulation is then apphed to the calculation of 
conductance (section |V|) and density of states (section [VI| ). We conclude in section |V11 . 



II. DEFINITION OF THE MODEL 

We begin this section by upgrading the above qualitative introduction of the sublattice 
system to a more quantitative formulation. Quantum transport through the bulk of the 
system is described by the Hamiltonian 

H = Y.i,{tiv5^'' + Riv)^i','^ (2) 

{W) 

where creates a spinless electron at site / in state /i = 1, . . . , iV, the sum extends over 
nearest neighbour sites, and Rui are N x N Gaussian distributed random hopping matrices 
with moments 

^2 

{RlvRw^ ) = J^^fiu^pi'u', A < 1. 

Apart from the Hermiticity condition Rut = R\,i matrices on different links are statistically 
independent. These random matrices compete with the regular contribution to the hop- 
ping matrix elements, tui. To be specific, we set t/// = 1 + (— )a if the smaller of the two 
neighbouring site indices / and /' is even (odd). The real parameter a is a measure for the 
staggering strength. Notice that tu = 0{1) ^ A implies that we are dealing with a weakly 
disordered system. 

At both ends, a number of sites is coupled to leads (see figure |^). Quantum propagation 
within these leads is assumed to be generic, i.e. governed by a Hamiltonian without sub- 
lattice symmetry. To describe the coupling, we add to our bulk Hamiltonian a tunnelling 
contribution 
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1=1,2,... 

H^= E (cXad^ + h.c. 



1=1,2,. 
1 

l=L,L-l,... 

where d^^^^^ creates an electron propagating in the left (right) lead in a certain state a = 
1, . . . ,M ^ 1 and we have introduced a composite index a = comprising site and 

orbital index of the bulk system. The coupling matrix elements W^^^ are subject to the 



orthogonality relation |29 



E KiL-i),wii_,)^^, = fmtSw, (3) 

where f{x) is some envelope function decaying on a scale of 0{1) and normalized through 
J2if{^) = 7 ^ 1- The function / and parameter 7 describe profile and strength of the 
coupling, respectively. Why did we introduce the multi-site coupling operators (H) instead 
of connecting just the two terminal sites I = 1, L to the lead continuum? Modelling the 
coupling in a more general way is motivated by the presence of the site number staggering 
phenomena mentioned in the introduction. The above implementation of the coupling op- 
erator is sufficiently flexible to selectively probe these effects (sections [V A| and [VI A| ) or to 



average over any boundary oscillatory structures (sections |V B| and |V1B| .) 

To conclude the definition of the problem, let us introduce the Green function, 

G{z) = iz- H + i7r{sgnlmz) ^ ly^iy^^j , (4) 

\ C=L,R J 

where W^W^'^ = {J2a^lli,aWa,ii'i} is an operator describing the escape of electrons from 
the bulk system into the leads |2^J5^. Expressed in terms of these objects, the Landauer 
conductance of the system assumes the form 

M 

g = i2nY Y: KaKbK'^a'a (G'.;3(0+)G^w(0-)) . (5) 

ab=l 

Our second quantity of interest, the density of states per site, is given by the standard 
expression 
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Ke) = -— Im5:(G.„(e+)). (6) 

Finally, to prepare the field theory formulation, let us consider the symmetries of the 
problem. Expressed in the notation introduced above, the chiral symmetry assumes the 
form (3"3]_|_ = 0, where (J3 = {{—YSu'} and H denotes the bulk part of the Hamiltonian. 
(Coupling to a non-sublattice continuum breaks chirality.) The presence of this symmetry 
implies invariance under the two parameter family of transformations 

cj cje~^^, Q e~^^ci, I even, 

cj c^e+^^ Cl e+'^ci, I odd, (7) 

where 2:1^2 are complex numbers. Depending on the choice of these parameters, (|^ expresses 
the standard U(l)-invariance of a model with conserved charge {zi = —Z2), or the axial 
symmetry characteristic for chiral systems {zi = ^2)- (For the time being we treat the trans- 
formation as purely formal i.e. ignore the fact that for a general choice zi^2 the transformed 
operators are no longer mutually adjoint.) We will come back to discussing the role of these 
symmetries after the effective field theory of the system has been introduced. 



III. FIELD THEORY 



The construction of the low energy effective field theory of the sublattice wire follows 



the standard route of deriving nonlinear cr-models of disordered fermion systems |31|, there 
are no conceptually new elements involved. Referring to Appendices ^ and ^ for technical 
details of this derivation, we here discuss structure and key features of the resulting model. 

As shown in the Appendices, the field theory representation of conductance and DoS is 
given by 

9 = - {{T - T-^),2(0)(T - T~%,iL)) , (8) 

and 
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^Re £dr(n,ir) + T,-,\r)), (9) 



respectively. Here, the bulk DoS, uq = T(r) is a field taking values in the supergroup 
GL(1|1) and the continuum variable r G [0, L] replaces the site index. The angular brackets 
stand for functional averaging 



{...) = J DTe-^[^l (. . .) 



over a functional integral with action S = Jq drC and effective Lagrangian 



C = C{l + Ctop + Ct + Cz + ^Ga^de- (10) 



The individual contributions are given by 

Aop = -^^^str (TdrT-'] 

£fl = -4str {drTdrT-'), 
Id 

= -^'-^ sir (T + T-^), 



= str (T(r) + T-\r)) [5{r) + 5{L - r)] 



StT{TdrT~^) , (11) 



where C is a coupling constant that need not be specified other than that it is small, 
C = (9(1) ^ (A^, M) and the parameter / = measures the degree of staggering. (Notice 
our / is identical to the control parameter / defined in Ref. |^.) Finally, we have introduced 
a parameter 

which will later identify as the localization length of the system. 

To prepare the further discussion of the functional expectation values, let us briefly 
discuss the internal structure of the field theory. We first note that save for the values of the 



coupling constants, the structure of the action (TU) can be anticipated from inspection of Eq. 
(0) and its supersymmetric extension: the field theory approach starts out from a promotion 
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of the fermionic operator representation to a supersymmetric formulation in terms of 
Bose and Fermi fields. Within the latter representation (cf. e.g. Eq. ( |A1| )), the space of 
permissible invariance transformations is enlarged. The formerly complex parameters e^' are 
replaced by two-dimensional matrices Tj acting on the bosonic and fermionic components 
of the theory. Any sensible subsequent manipulation done on the functional integral must 
respect this invariance property. On the level of the effective theory described by S[T], the 
transformation acts as T — T1TT2 and indeed one verifies that the contributions £top, -^a and 
C(f) of ([lO|) are invariant under this operation. Further, the two building blocks str {TdT~^) 
and str ((9T(9T^^) are the only operators with < 2 derivatives compatible with the global 
GL(1|1) symmetry. In other words, the gross structure of the bulk action readily follows 
from the invariance criterion. For finite energies or coupling to the leads chirality is broken 
and global U(l) remains the only symmetry of the Hamiltonian. Within the supersymmetry 
formulation, the set of allowed transformations is then reduced down to configurations with 
Ti = (the super-generalization of U(l)). The operator str (T + T~^) is the minimal 
positive (see below) choice compatible with the restricted symmetry. Summarizing, the three 
terms str {TdT"^), str {dTdT~^), and str (T + T~^) exhaust the set of operators compatible 
with the global transformation behaviour of the model. 

As for the coupling constants of the theory - not specified by symmetry arguments but 
all derived in Appendix ^ - notice that contrary to the standard a-models of systems with 
WD symmetry two, instead of just one second derivative operators appear in the action. 
Mathematically, the reason for the appearance of two contributions is that the target mani- 
fold of the theory, GL(1|1) is a reducible symmetric space; it decomposes into two irreducible 
factors, a point discussed in detail in Ref . . Each of these factors can be endowed with its 



own metric which implies the existence of two independent second derivative operators in the 
theory. Physically, the presence of the non-standard operator ~ str^(T(9T~^) has profound 
consequences for the behaviour of the 2(i- version of the field theory 0: the RG analysis of 
the model shows that the coupling constant of this operator grows under renormalization 
while driving the coupling of the energy operator ~ str (T + T~^) to large values. At the 
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same time the coupling constant of the standard gradient operator ~ str {dTdT~^), essen- 
tially the conductance, remains un-renormalized. In one dimension, the situation is different. 
Contrary to what one might expect, the contribution iZcade is not remotely as important as 
in two dimensions, and it is another operator that drives the localization behaviour of the 
model. In parentheses we remark that the target space of the transfer matrix approach to 
the problem, GL(M)/U(M) ~ SL(M)/SU(M) x 7^+ factors into two components, too 



Accordingly, the Fokker-Planck equation governing the 'Brownian motion' on that space 
is controlled by two independent coupling constants, both determined by the microscopic 
definition of the model. The second of these contributions, essentially the analogue of our 
Lagrangian C^, leads to non-universal corrections to the overall picture which perish in the 
limit of a large number of channels. 

A second aspect discriminating the present model from its WD relatives is the appearance 
of a first order gradient operator in the action. In fact, the presence of this contribution 
might raise suspicion: although allowed by symmetry, str {TdrT~^) is not invariant under 
space reflection r — r, in contrast to the microscopic parent model (for a = 0.) The 
resolution of this puzzle lies in the fact that Stop = / £ is of topologcial origin and, although 
not manifestly so, does respect the space inversion property. We will discuss this point 
momentarily after the internal structure of the field manifold and the integration measure 
have been specified. 

Functional integrals can only be defined on manifolds that are Riemannian, i.e. endowed 
with a positive metric. The supergroup GL(1|1) (like all other supersymmetric spaces that 
appear in the context of field theories of disordered Fermi systems) does not fulfil this 



criterion, a point discussed in detail in Ref. [|T8|. However, it does contain a maximally Rie- 



mannian sufemanifold Ai, defined as follows: the boson-boson block of A^, a one-dimensional 
manifold by itself, is isomorphic to the non-compact symmetric space GL(1)/U(1) ^ TZ^, i.e. 
the positive real numbers. (This space is trivially Riemannian.) The fermion-fermion block 
is isomorphic to the compact symmetric space U(l). No limitations in the Grassmann valued 
boson-fermion sectors of the theory are needed since the whole issue of convergence does not 
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arise here. Summarizing, A4 = TZ'^ x U(l), where the notation is symbohc, specifying the 
boson-boson and fermion-fermion sector, respectively. 

A convenient field representation respecting these convergence criteria is given by 

T = kak~^, k = exp | ] , a = exp | ] , x,y eTZ. (12) 

V'^ / \ wJ 

In this parameterization, the group integration / dT extends over the degrees of freedom 
X, y, 7], u, without further constraints. The invariant group measure associated to the repre- 
sentation (|12D are defined in Eqs. (|ClD and ( p2D . 

We are now in a position to discuss the role of the contribution Stop- First notice that for 
sufficiently strong lead coupling the boundary action projects the fields onto the group 
origin T{x,y) = U, i.e. enforces x(0) = x{L) = and y{0) = 27ck, y{L) = 2TTk', where k, k' 
are integer. As a consequence, the first derivative operator can be written as 

N + f 



StoAT] = Y~ JdrstT{TdrT 

N + f 



2 

N + f 



drdr str In T 



str (In T(L) - In T(0)) 



2 

= i7r{N + f){k-k'). (13) 

This makes the topological nature of the term manifest: it counts winding numbers in the 
fermionic sector of the theory. The integer k — k' is a. topological invariant characterizing 
each field configuration T(r). Further notice that for / = 0, 

g-5top[T] ^ ^+iTTN{k-k') ^ ^-inN(k-k') 

Since it is the exponentiated action and not the action itself that matters, the last identity 
tells us that the global sign of the first gradient operator is irrelevant (all for / = 0). This 
settles the above raised issue of the behaviour of the model under space refiection: although 
the first order derivative is not invariant under r — — r the exponentiated action is. 
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IV. TRANSFER MATRIX APPROACH 



Principal aspects of the system properties we are going to analyse are non-perturbative, 
i.e. cannot be obtained as power series in the coupling constants of the action. Progress 
with such type of problems can be made by applying the cr-model transfer-matrix technique 
an approach conceptually similar to the DMPK formalism. 
We begin by defining the two functions 



T(r)=T2 



T(L)=T2 

Expressed in terms of these objects, the DoS assumes the form 

^(e) = ^ Re dr j dT ^^(1, T, r) (Tn + Tf/) F^(T, U, L - r), (14) 

where we have used that for sufficiently strong coupling to the leads, the boundary config- 
urations T are close to unity. (Throughout much of this paper we will consider the DoS 
of coupled systems. For sufficiently large systems, the choice of boundary conditions is 
inessential, a point to be verified below.) Similarly, the conductance obtains as 



9 



^g^^str(T+T-)(y - r-l)i2 Yl{T,T',L) {V - r-%,e~^^'^^^'+^"'\ (15) 



From Eqs. (|T4D and (|T5D it is clear that the functions Yl^r encode the essential system 
properties we are interested in. 

As a first step towards the computation of these objects let us explore how the sym- 
metries of the action translate to symmetries of Yl ji. We first consider the case z = 0, 
relevant to the analysis of the conductance, where the action is fully invariant under GL(1|1)- 
transformations. The invariance C[T] = C[TiTT2], T12 = const, then directly implies 
Ylji{T,T' ,r) = YL ji{(TiTT2), (TiT'T2),r). From this identity one readily verifies that 
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Yr[T, T', r] = Yr[TT'-\ U, r] = Yn[T'-'T, U, r], 

YL[T,T',r] = YL[tXT-\r] = YL[l,T-'T',r], {z = 0). 

In other words, for z = the arguments of the heat kernels enter in an invariant product 
type form and it suffices to consider the reduced functions 

r^(T,r) = Yr{T, U,r), r,.(T,r) = 1^(1, T,r), (16) 

depending on a single argument only. The same invariance property (now evaluated for 
T2 = Ti^) implies that YiniT^r) = YijiiTiTTi"^ ,r). Imagining the argument T to be 
represented in the polar decomposition (|T2|) and setting Ti = the argument can be 
reduced to the diagonal matrix of eigenvalues a: 

YL,R{kcik-\r)=YL^R{a,r). (17) 

For -2 7^ 0, the invariance of the theory collapses down to transformations with T2 = T{^. 
However, the representation of the DoS above implies that we are solely interested in func- 
tions of the type (|16|), with second argument set to unity, anyway. Since transformations 
T TiTTi^ are still permissible, these objects depend on the eigenvalues of the argument 
matrix only, as before for the case z = 0. Summarizing, in the analysis of both conductance 
and DoS, it is sufficient to consider functions Yl,r depending on a single argument with 
invariance property (0). 

Following the basic philosophy of the transfer matrix approach, we will compute the 
functions Yl^r{T, r) iteratively, by asking how much they change under infinitesimal variation 
of the arguments r — r + e. Considering the function Yl for definiteness, we first notice 
that, by definition, 

F^(T,r + e) = J DTe-/r^'^^'^[^V,.(T(r), r). 
For sufficiently small e, the action can be expanded and we obtain 

^^iv+/str(rT-)+^^str(T+T-) (y^(r,r) - YL{T,r)) , 
20 



where we have used that, due to supersymmetry, J VT exp{— J^'^" C[T]) x 1 = 1. We have 
also temporarily set the coupling constant of the operator Ccade to zero. As mentioned 
above, this term does not have much relevance in the present context. We will briefly 



discuss its role in section VE 



From hereon, the derivation of an evolution equation for Yl is conceptually straightfor- 
ward: the exponential weight factor in the first line of the equation enforces that T' is close 
to T, symbolically, T'T^^ = 11 + 0{e). We should thus expand T' around T and evaluate 
the integral / dT' perturbatively in e. This expansion is most conveniently done in the po- 
lar coordinates introduced above (because the heat kernel depends on the radial degrees of 
freedom, only.) As a result of a calculation similar but much more simple than the one for 



the standard a-models with their larger matrix fields ||31| one obtains the Schrodinger type 
equation 

[Tdt - 4(D ± A)2 + y)) Yl^ (a, r) = 0, (18) 
where V{x, y) = — is (cosh (x) —cos (y)) and we have introduced the dimensionless parameters 

s = TTU^e. 



r 



Physically, t is the length coordinate measured in units of the localization length ^ and s the 
energy measured in units of the single particle level spacing = of a system of length ^. 
Further, the symbol D = {D^, Dy)^ denotes a vector differential operator defined through 



where the constant vector 



Dx = dx coth (- — 

2 V 2 

Dy = dy+'^- coth (^L_pL 



(19) 



(20) 



Finally, evaluation of (18) for asymptotically small times t \ leads to the initial condition 



\imYL,R{x,y,t) = S{x,y) = lime 8t(^^+s^^). 



(21) 



21 



It is very instructive to interpret the structure of the evolution equation ([TSD in the 
hght of the analogy between field theory and point particle quantum mechanics on GL(1|1) 
discussed in the introduction. Within this picture, the functions Yl[T, r] acquire the meaning 
of Green functions, i.e. transition amplitudes for the propagation between the origin of the 
manifold T = U at time and a final configuration T at time t. The evolution equations 
(|T8|) becomes a a time-dependent Schrodinger equation with quantum Hamiltonian, H = 
-2(D - A)2 + V{x,y). While the term V{x,y) = -is(cosh(x) - cos(y)) = -f str (T + 
T~^) simply represents the potential inherited from the Lagrangian, the 'kinetic' operator 
is more interesting. The covariant structure (D — A)^ describes minimal coupling to the 
constant vector potential where the unfamiliarly looking structure of the derivative operator 
D is a consequence of the fact that our particle lives on a curved manifold. Indeed, it is 
straightforward to verify that for A = 

D ■ D = 5] J-'d,Jd,, 

i=x,y 

where J{x, y), specified in Eq. ( \C2\j , is the square root of the determinant of the metric tensor 
on GL(1|1). This structure identifies D ■ D as the radial part of the Laplace operator on 
GL(1| 1). ('Radial', because the two coordinates x and y, spanning a maximally commutative 
sub- algebra of the Lie algebra of GL(1|1).) 

Summarizing, we have identified H as the Hamiltonian of a charged particle on the group 
manifold GL(1|1). Our next task will be to compute its Green functions y^/ij- We begin 
by considering the case of a free particle, = e = 0. The solution of this problem will 
contain the information needed to compute the conductance. 

V. CONDUCTANCE 

In this section the general formalism developed above is employed to compute the con- 
ductance g of the system. Our main objective will be to understand the impact of topology 
on the localization behaviour of the system. However, before embarking on this analysis 
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it is tempting to digress for a moment and to briefly consider the transport behaviour of 
short systems, specifically the aforementioned anomalous sensitivity to the coupling to the 
leads. Being not directly related to the mainstream of the paper, the technicalities of this 
discussion have been deferred to Appendix |^ and we here restrict ourselves to a summary 
of results. 

A. Digression: Conductance of Short Wires 

In the present context, the phrase 'short' means that systems in the quantum dot regime 
L < ^/{M'-f) are considered: the conductance is not so much determined by the bulk 
transport properties of the system but rather by the strength of the coupling to the leads. 
Moreover, and this is a speciality of the sublattice system, the parity of the coupling, i.e. 
the even- or oddness of the connecting sites, turns out to play a crucial role. More precisely 
we find that (a) for systems where at both ends a number of sites of alternating parity are 
connected - the setup considered in much of this paper - the short system conductance is 
given by 

9 = ^, 

in accord with the behaviour of conventional quantum dots. The same result obtains (b) 
for systems where only one site at either end is coupled and these sites have opposite parity 
(even/odd or odd/even). However, (c) for single site coupling with even/even or odd/odd 
connectors, the conductance vanishes. To heuristically understand the phenomenon, it is 
instructive to consider the toy-model case of a strictly one-dimensional clean sublattice 
system. For zero energy, the wave length of current carrying excitations is commensurable 
with the lattice spacing. This means that the relevant quantum wave functions have nodes 
at alternating sites. E.g., for an even-even configuration a state entering from the right 
has zero quantum amplitude at the exit site on the right. This implies a total blockade of 
electric current. A less obvious fact is that this phenomenon survives generalization to many 
channels and disorder. 
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We repeat that all these results are obtained for short systems; field fiuctuations, de- 
scribing the propagation of spatially non-uniform diffusive excitations, are neglected. An 
interesting question, not considered in this paper, is how such modes would affect the con- 
ductance as the system size is increased. 

B. Reduction of the Problem to (0-dimensional) quantum mechanics on GL(1|1) 

We next turn back to the discussion of large systems (equihbrated coupling of type 
(a) understood.) Inspection of the basic expression ([T5|) shows that the problem factorizes 
into doing the boundary integrals and analysing the bulk behaviour of the heat kernel, 
respectively. 

We begin by discussing the boundary regions. Following a line of arguments developed 



in Ref. p2[, we first notice that due to the exponential weights ~ exp(— str (T + T ^)), 



the integrands are confined to the immediate vicinity of the group origin. This suggests to 

write T = exp W, and do the integrals over generators in a Gaussian approximation. 
fu a\ 

Setting = and expanding in coefficients we arrive at integrals of the structure, 

g = const. X J dudvdpdae-^'''''+'''+^''^^aF{T{u,v,a, p)), 

where the symbol F{T), represents the functional dependence of the heat kernel on the 
boundary field. Evaluation of the Gaussian superintegral leads to 

g = const. X :^^dpF{T)\^^-^. 

Doing the same procedure for the second boundary integral and fixing factors we obtain 

9 = dpd^,YL{T,T',L)\^^^,^-ii_ = dpdMTT'-\L)\^^T,^-^, 



where the second expression contains the one-argument heat kernel introduced in Eq. (p!6[). 
According to this expression, the conductance is obtained by second order expansion of the 
heat kernel around the origin. We next note that due to the invariance property (JT^ the 
expansion starts as 
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F^(f = e^,L) = l + cistT{W) + C2stT{W^) + C3[stT{W)f + .... (22) 
Now, in our case, 

/ / 0\ / -a'W (\a'p -a' 
W = \n{TT'-^) = In exp exp = 

V Vp / Vo J) \ P \<y'P, 

Substitution of this expression into ( p2D and differentiation leads to 

9 = 2C2, (23) 

i.e. the problem has been reduced to fixing the coefficients of the series expansion of Y^. 
Notice that this series representation is totally determined by the features of the bulk system; 
all aspects related to the coupling to the leads have disappeared from the problem. 

We now have to face up to the principal task, the calculation of the heat kernel. 
Its interpretation as the Green function of the problem suggests to begin by represent- 
ing Yl through a formal spectral decomposition: consider Eq. ([T8| ) at zero potential, 
(— 9j — 4(D — A)^) Fi(a;, y, t) = 0, and suppose we had managed to find a set of eigen- 
functions, 

-4(D-A)2vl/W = e„vl/M. (24) 
Assuming completeness, we can then span the heat kernel as 



where the expansion coefficients are determined through the initial condition (|2lD . Provided 
the spectrum is suitably structured (positive and gapped against some low lying levels), and 
keeping in mind that we are interested in asymptotically large values of t, the sum may be 
restricted to a limited set of n's. The problem thus reduces to (a) exploring the low-lying 
spectral content of the operator D — A, and (b) fixing the expansion coefficients. 



25 



C. Spectrum and Eigenfunctions of the 'Kinetic Energy Operator' on GL(1|1) 



In spite of the unfamiliarly looking coordinate representation of (D — A)^ analytic 
progress with the problem is straightforward, the reason being that this operator is noth- 
ing but a plane wave operator in disguise. To make the hidden simplicity of the problem 
manifest, let us first remove the dependence of the differential operator on the (pure gauge) 
potential A: transformation \E'„(x,?/) —>■ e^^''^'~*^''\E'„(x, = ^„(a;,?/) brings the eigenvalue 
equation (P^ ) into the form 

where A = DD is the radial part of the Laplacian on GL(1|1). Notice that the structure of 
this equation does not imply that the vector potential has disappeared from the problem. 
It has merely been transferred from the differential operator to the boundary conditions 
attached to the differential equation. While irrelevant in the non-compact sector of the 
theory, changes of the boundary conditions in the compact sector generally cause qualitative 
effects. To appreciate this point, notice that the un-gauged Hilbert space of the problem has 
periodic boundary conditions, "^nix, y) = '^n{x,y + 2'K). Yet, after the gauge transformation, 
= ^n{x,y + 27r)(-)^+^, i.e. for odd or non-zero staggering a transmutation 
to twisted boundary conditions has taken place. Needless to say that this change bears 
consequences for the spectral structure of the problem and, therefore, for the transport 
behaviour of the system. 

To make further progress, we subject the eigenvalue problem to the similarity transfor- 
mation 

A^j'/'^AJ-'/^ = A = dl + dl 
This change of representation entails an enormous simplification of the problem. The un- 



invitingly looking operator A has become a fiat two-dimensional Laplace operator [33 



However, notice that the transformation by J^^^ effects yet another change in the bound- 
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ary conditions: due to J^/^(a;,0) = — J^/^(x, 27r), the transformed states are subject to the 
condition 

^„(x,y) = §„(x,y + 27r)(-)Wi. 

At this point, the solution of the eigenvalue problem has become a triviality. The equation 

-4{dl + dl)^n{x, y) = en^n{x, y), 

is solved by the exponentials ^k,i{x,y) = e*^'^"'"*^*^, where {k,l) = n are two 'quantum 
numbers', pk^i the associated momenta and the eigenvalues eki = 4(p| + pf). From these 
states we obtain our un-gauged and un-transformed original wave functions as 

^,,(^,1/) = sinh (^^) e(*^-'^)^+^(^^+'^)^. (25) 

To give this set of solutions some meaning, we need to specify the range of permissible fc's and 
Vs. In the compact sector the situation is clear - the circular boundary condition specified 
above enforces pk = k — ^y^, with halfrntegei k. The conditions to be imposed in the 
non-compact sector are tightly linked to the integrability properties of our wave functions: 
The space of radial functions on GL(1|1) is endowed with a natural scalar product, viz. 

/oo /•27r 
dx / dyJ{x,y)f{x,y)g{x,y). (26) 
-oo Jo 

We demand that the eigenfunctions contributing to the spectral decomposition of the heat 
kernel be square integrable w.r.t. ( , ). Inspection of Eq. ( p5D shows that this requirement 
enforces pi = I — i^^-^, where I may be arbitrary and real. 

Finally, we need to specify a set of functions sufficiently complete to generate an expan- 
sion of the heat kernel. The present problem does not come with a natural Hermitian or 
symmetric structure (i.e. for finite A the kinetic energy neither symmetric nor Hermitian.) 
However, defining a 'fake complex conjugation' through ^'^^ = it is straightforward 

to show that 

{^ki,^k'i') = {271^6(1 -l')6kk'. (27) 
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We may thus attempt to represent sufficiently well behaved (for the cautious formulation, 
see below) functions as 

g{x,y) = E / dl gki^ki{x,y), gui = {2n)-^{^M,g). (28) 
k 

Before applying this procedure to the heat kernel, let us summarize our main findings for 
clarity: the radial Laplacian on GL(1|1) is diagonalized by the set of functions 



^ki{x,y) = sinh ( ] e'^'^^'y^ . (29) 



where k E Z + 1/2, I real and the eigenvalues are given by 



eki = 4 



(30) 



(Notice that the appearance of an imaginary part proportional to the strength of the vector 
potential is due to the fact that for finite A, the Hamiltonian of the theory is neither 
symmetric nor Hermitian.) The spectral decomposition of radial functions is defined through 
Eqs. (g6D, (grp and (|8|). 



D. Computation of the Conductance 

We now apply the machinery developed in the last section to the analysis of the heat 
kernel. In principle, the strategy seems to be prescribed by what was said above. We should 



determine the Fourier coefficients of the initial configuration (0), 6ki = (^jfe, 5), from where 
the heat kernel would follow as YL{x,y,t) = {2n)^'^J2kI dl 5kie~'^'^^^'^ki{,x,y). There is a 
problem with this procedure, viz. the expansion coefficients of the '(5-distribution' 5(x, y) 
vanish. Indeed one verifies that in the limit t — 0, the support of the Gaussian in (pT]) shrinks 
to zero while its maximum remains limited by one. This readily implies (^E', 5) = 0. The 
reason for this pathological behaviour is that our radial theory memorizes that it derived from 
a supersymmetric parent theory. In a sense, supersymmetry can be interpreted as a theory 
on a zero- dimensional background, i.e. there is no singular 'volume factor' compensating 
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for the vanishing support of the (5-function, as would be the case in spaces with positive 
dimension. 



The problem can be circumvented by a cute trick Instead of Fourier expanding Y^,, 
we consider the function Yl — 1. Since unity by itself solves the heat equation, no harm 
has been done and all that has changed is the boundary condition: lim(^o(^L(3;, y, t) = 
S{x,y) — 1, a function that equals minus unity almost everywhere save for the origin where 
it vanishes. 

We thus represent the heat kernel as 

Yl{x, r) = 1 - (27r)-2 J2 [ dl Ui^kiix, y)e-'-\ (31) 

k 

where 



i{lx+ky) 



/ — ik 

(To obtain the last equality, it is convenient to first do the x-integral. Closure of the 
integration contour in the upper/lower complex half plane for positive/negative / yields a 
semi-infinite sum over residues of the sinh-f unction, along with a y-integral that is of simple 
plane wave type. Doing sum and integral one obtains the result.) 

Substitution of this result into the expansion of the heat kernel now yields 



Yl{x, y,t)-l = ^ dl sinh 

vr ^ J L — ik 



k 

where in the last line we have switched back to a coordinate invariant representation. Com- 
parison with Eqs. (|2^) and (p3D finally leads to the identification 



1 coo 4,L 

^" kez+i/2 



oo 
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where we have inserted the exphcit form of the eigenvalues. Notice that all manipulations 
leading from the original cr-model representation to the above Gaussian integral representa- 
tion have been exact. 

We next evaluate this result in the two limiting cases of physical interest, L <^ ^ (Ohmic 
regime) and L ^ ^ (localized regime.) Beginning with the Ohmic case, we first notice 
that for L <C ^ many terms contribute to the k summation implying that the sum can be 
approximated by an integral. Thus, 

16L 



^ 2 V47rL ; J 161 



As for any ordinary conductor, g is inversely proportional to the system size; the parity of 
the channel number does not play a role. 

In the opposite case, L ^ ^, only those discrete indices that minimize the exponent 
contribute to the sum. Specifically, for an even channel number and no staggering, 

9 ~ I ^ J ^ ^ ' ^ ^ °" 

In contrast, for odd and / still zero, the exponent vanishes for the half integer k — ^ 

UAV'^ A^odd,/ = 
^ 2 \AttL y ' 



depends algebraically on the system size. Finally, it is clear that for non- vanishing staggering, 
/ 7^ 0, intermediate types of behaviour are realized. E.g., an N even chain with staggering 
/ = ±1 behaves like a symmetric odd chain, etc. 

E. The Role of the Gade Term 

Before leaving this section let us briefly discuss the role of the, so far neglected, Gade 
operator S'Gadc[7']. The inclusion of this term in the derivation of the heat equation is 
straightforward. As a result, the planar Laplacian A — + dy gets replaced by 
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where ^ ^ 1. This operator is still diagonalized by the plane waves discussed 



above. The eigenvalues change to 
eki =4 



Recapitulating the computation of the conductance, one finds that the small dilatation 
introduced by finite values of fj does not affect the long range transport behaviour of the 
system. 

All this is compatible with the structure of the DMPK transfer matrix approach to the 
problem. As mentioned above, the DMPK evolution equation is controlled by two coupling 
constants. One of these, in Ref. |TD[ denoted by 77, is small, 77 ~ 0{M~^) and becomes 



inessential in the limit of a large number of channels. The analysis above suggests that 
the coupling constant of the Gade operator, oc r/, and the r] of the DMPK approach play 
the same role. This analogy is supported by the above discussed geometric structure of the 
target spaces of the two theories. 

VI. DENSITY OF STATES 

We now turn our attention to the low energy density of states of the sublattice system. 
As with the conductance, we will first consider the behaviour of short wires, and then discuss 
the localized regime. 

A. Density of states of short wires 

As in section |V A| we consider a short sublattice wire of length L < ^/{M'-f). The 
spectrum of such systems exhibits structure on the scale of the mean level spacing. In 
order not to blur these fine structures, the coupling to the external leads will be switched 
off throughout this section. Following the logics of section |V A| , one would then conclude 



that the action reduces to S[T] = J drC.^[T] = — 'j|str(T + T ^), where the matrix T 
parameterizes a zero mode configuration, s = ^ and A = (z/L)~^ is the level spacing of the 
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isolated system. (Temporarily deviating from the convention of the rest of the paper, in this 
section s measures the energy in units of the total level spacing and not the level spacing of 
a localization volume.) This presumption is almost but not quite correct. The point is that 
our so far discussion of the low energy action implicitly assumed that the number of sites of 
the system is even. In the opposite case, an extra contribution S'top,2[^] = A^str ln(T(L)), 
derived and discussed in section |B 1| , appears. The structure of this term reflects the fact 
that due to the presence of one uncompensated site the global GL(1| l)-invariance of the 
model is lost. {Stop,2[T] is not invariant under T T^TTr even for constant Tl,r-) In the 
majority of cases this extra term is of little interest. However for the spectral properties of a 
short isolated system, the presence of S'top,2 bears crucial effects. Indeed we will see that this 
term is responsible for the formation of staggering phenomena akin, and probably related 



to the effects discussed earlier in section |V A 



Adding the two contributions to the action and using Eq. (|) we obtain 

for the cr-model representation of the zero- dimensional DoS. The task thus is to integrate 
over a single copy of the target manifold GL(1|1). For sufficiently small energies e ~ A, 
the action is not large enough to confine the integrand to the origin of the group manifold 
implying that the integral has to be done non-perturbatively. Referring to Ref. and 
Appendix ^ for technical details we here merely display the final result of this integration 
procedure. 



i^[S) = —r 
^ ' 2A 



J^{s) + Jt{s) 



L even 



u{s) = N6{e) + — [jUs) - Jn+i{s)Jn-i{s)\ , L odd (32) 

The structure of these DoS profiles is shown in Fig. ^ for the example = 3. 

Eqs. (^) have been obtained earlier within pure random matrix theory |T^, and its 



supersymmetry implementation [P^,pT|]. Studies of chiral random matrix ensembles were 



largely motivated by the relevance of the spectral structure of effectively zero dimensional 
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chiral systems in finite size lattice QCD. More generally, 'microgaps' of the type shown in 
Fig. are an omnipresent side effect seen in the spectrum of generic chiral systems with 
finite mean level spacing A. 




s 



-6-4-2 2 4 6 

FIG. 3. Density of states of the short sublattice system. Solid line: DoS of a system with L 
even. Dashed line: DoS of a system with L odd and = 3. (The singular zero energy spike 
causing the spectral depletion up to values s ~ 3 is not shown.) 



Qualitatively, origin and structure of these gaps can easily be understood. First, the 
chiral symmetry crs]^ = entails that for any non-vanishing energy level e„ its negative 
— e„ is an eigenvalue, too. Disorder generated level repulsion prevents these states from 
coming close to each other (on the scale of the mean level spacing) which explains the 
presence of the spectral gap in the L-even case. For L odd, this picture has to be modified. 
To understand what is happening, let us imagine the Hamiltonian as a block off-diagonal 
matrix in sublattice space: 

Z\ odd 

Since L is odd, the number of odd sites exceeds the number of even sides by one, i.e. the 
blocks Z are rectangular with N{L — 1)/2 rows and N{L+ l)/2 columns. Now, any block-off 
diagonal matrix with k rows and / columns has |fc — / | eigenvalues 0. Applied to our system, 
this means that the L odd system has zero modes for any realization of disorder (the 



H 



even. 



5-function contribution in (|32D .) Other levels repel from this concentrated accumulation 
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of spectral weight which explains the bathtub type suppression of the DoS up to energies 

Summarizing we have seen that the L even/odd staggering behaviour observed earlier in 
connection with the conductance pertains to spectral properties. Above we had argued that 
the vanishing of the conductance in the L even case was due to the peculiar spatial profile 
of zero energy wave functions. In view of (|32D it is tempting to relate the same effect to the 
vanishing of spectral weight at zero energy, although we have not analysed this picture any 
further. 

Keeping in mind the tendency towards buildup of zero energy spectral weight in systems 
with mis-matched sublattice structure we next turn back to the analysis of large systems. 



B. Heat kernel for finite energies 

The aim of this section is to compute the DoS of a large system with L ^ ^. To 
facilitate comparison with the behaviour of the conductance we will drop the assumption of 
isolatedness and again couple the system to leads. Notice that this stands complementary to 
the analysis of Ref. where a closed system was discussed. That we will obtain identical 



scaling of the DoS is proof of the (in view of the existence of topological zero modes not 
entirely obvious) assertion that boundary conditions do not affect the bulk spectrum. 

We will not be able to compute the DoS z/(e) for arbitrary e. Instead, an asymptotic 
expression valid for low energies will be derived. Remembering the behaviour found for 
small systems we anticipate nonanalytic behaviour, u{e) ~ ln"(|e|)|e|, where the logarithmic 
factor is due to potentially existent localization corrections to the zero mode behaviour. Our 
objective is to identify the most singular contribution of this type. 

Starting point of the analysis is the transfer matrix representation ( p^ of z^(e). To 
evaluate this expression we need to compute the functions Y^^^, now for finite potential 
V. We will do this following a procedure developed in Ref. Starting out from a 



spectral decomposition of the type (^Tj), we first notice that only zero energy eigenfunctions 
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contribute to long distance behaviour of lL/_R(r). This means that our 'time' dependent 
'Schrodinger equation' (|18|) can effectively be replaced by the stationary form 



t > e : (-(D ± A)2 + 7/(cosh(x) - cos(y))) Yl/r{x, y) = 0, (33) 

where we have substituted the explicit form of the potential and omitted the spatial argument 
of 1l/_r for simplicity. For later convenience we have also analytically continued from real 
to imaginary energy arguments, —is'^ ^ At] > 0. Substitution of these functions into Eq. 
(|1^) yields the reduced representation 

u{E) = ^0 (|l + Re ^ / dxdy '-^^^^^^^M^, y)Ynix, j . (34) 

To obtain this equation we have expanded the pre-exponential sources in Grassmann vari- 
ables and integrated over these. The constant contribution uq appears as a consequence 
of the Efetov-Wegner theorem. (The representation of the pre-exponential term in polar 
coordinates contains a purely non-Grassmann contribution. Integration over a term of this 



type obtains the integrand at the origin [31] which, in our case, equals unity.) 

We now turn to the actual computation of the functions Y^/ji. The first step is the 
derivation of a set of matching, or boundary conditions relating the heat kernel to its ?7 — >■ 
asymptotics. To this end we evaluate the t] = spectral decomposition (^) in the limit 
t ^ ^. Neglecting contributions that decay exponentially in t/^, it is straightforward to 
obtain the asymptotic expressions: 

Yl (x, y) — > 1, even, 
y.(x,y)^i(e±(^-^^) + l), AT odd, 

where we have temporarily neglected the staggering parameter /. The structure of these 
functions will be motivated shortly. Turning to the case ?7 7^ 0, we next gauge and transform 
Eq. ( P5D as in section [V C| . Transformation YL{x,y) J^^'^{x,y)e^ 2^'^^"*^^yL (x, y) = 

It R 

Yl {x,y), where Vn = (1 + (~)^)/2 brings the equation into the form 

R 

{dl + dl - 77(cosh(a;) - cos(y))) fz./^(a;, y) = 0, (35) 
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while the transformed rj ^ asymptotics read 

YL/R{x,y) ^ sinh-^ (2^) 2 sgn (x)e-^snx(^)^ ^ ^^^^^ 
YL/R{x,y)'-^ coth(2^) l^'sgn(x), iV odd. (36) 

Here we have anticipated that the dominant contribution to the above double integral rep- 
resentation of the DoS will come from large values of the non-compact variable x. Eqs. ( pSf ) 
and ( PBD have the nice feature of full separability. Writing Y{x, y) = J\fYi{x)Y2{y) , where M 
is a normalization factor and the subscript L/R has been dropped for simplicity, Eq. ( pSj) 
can be traded for the set of decoupled ordinary differential equations 

{^l-le^-'^)Y^{x) = + O{v), 

The second equation is trivially solved by Y2{y) ~ e±«^ivf ^ xhe first line is a Bessel equa- 
tion. Its two solutions are given by /p^((2r7)-'^/^e'^'/^) and Kpj^{{2r]y^'^e^^^^'^). Discarding the 
exponentially divergent solutions Ii^ and using that for small arguments, Kq{z) ^ — ln{z) 
and Ki{z) ^ z~'^ + f (f ) 5 the normalization constants A/" can now be fixed by matching to 
the zero energy asymptotics in the limit of small r]. This obtains the approximate solution 

YL/Rix,y) '""l^^^ {8r]y/^Ki (^{27])^/^ e^^^/'^) e^^^nxf^ ^^^^^ 
YlMx, y) '^i'' -i-^K, ((2r/) V2eN/2) , N odd. 

It is now a straightforward matter to substitute this expression back into the above integral 
representation for v and to integrate over coordinates. The y-integration, extending over a 
purely harmonic integrand, is trivially done. (Notice that the sinh~^-factor in ( |3^ cancels 
against the factor (J^/^)^ from the similarity transformation). As for the x-integration, we 
note that due to the exponentially decaying asymptotics K^^z) ~ exp(— 2;/2) for \z\ ^ 1, 
the integral can be cut off at (217)^/2^1x1/2 _ 1 ^ | 

a; I ~ — ln(2?7). Within the domain of 
integration, the Bessel functions can be replaced by the small argument asymptotics specified 
above. Substituting these expressions it is then straightforward to obtain ij{ri) ^ z/q Re?] In^ 77 
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(even A^) and u{r]) ^ z/q Re (17 In^ r/) ^ (odd A^) for the small 77-asymptotics of the DoS. 
Analytic continuation back to real energies finally leads to the result 

z/(s) — z/o|s| In A^even, 



u{s) ^ - , , , , , AT odd, (37) 
s m s 



for the low energy behaviour of the DoS. Eqs. (37) agree with the results found earlier in 



Ref. m 



We finally discuss the extension of the above results to non-zero staggering. For non- 
vanishing / and even A^, the large argument asymptotics (|3^) generalize to 

YL/Rix,y) "^-^^^ 2e^(±^+^)(^-^^). (38) 

One can now follow the same steps as in the non-staggered cases above to obtain the result 

For non-zero /, the DoS vanishes in a more singular manner as in the non-staggered case. 
This behaviour provokes the question, how matching with the diverging profile in the case 
A^ odd, / = might be obtained. The principal structure of the theory entails that (A^ 
even// = 1) should be equivalent to (A^ odd// = 0). On the other hand, the / /" 1 version 
of Eq. (|39|) certainly does not agree with the divergent result (^Tf). 

To resolve this paradox, it is helpful to re-interpret the asymptotic expressions (|36|) and 
(^) within the quantum mechanical picture of the theory. Focusing on the compact sector 
and temporarily ignoring the factor J^/^ from the transformation to a flat Laplacian, these 
functions acquire the meaning of ground state wave functions \I'o of a one dimensional ring 
subject to a gauge flux For A^ even and zero /, an integer number of 'flux quanta' pierce 
the ring, and the ground state wave function carries zero 'momentum', ^^0(2/) oc 1 which, 
after multiplication with J^^^ leads to the first line of Eq. (p6|). In contrast, for A^ odd and / 
still zero, a half integer flux pierces the ring. This is a special situation in the sense that the 
ground state wave function is two-fold degenerate, i.e. '^o{y) = c+e*^/^+c_e~*^/^. Our earlier 

37 



analysis has fixed the a priori un-determined constants c± to a symmetric configuration. 
Inclusion of the non-compact variable and multiplication with J^/^ then leads to the second 
line of (H). 

We can now understand what happens as / is turned on for an even configuration: a 
fiux //2 is sent through the ring and the ground state wave function remains unique (cf. Eq. 
(^)) until f comes close to the critical value 1. In the immediate vicinity of the degeneracy 
point, the fact must no longer be neglected that our one-dimensional system is subject to 
a weak potential T]cos{y). For values of / such that the level splitting ~ 1 — / between 
the two nearly degenerate levels becomes comparable with the characteristic strength of the 
potential ~ r], the 'true' ground state configuration is given by the symmetric superposition 
of the two levels, as in the odd case. For these values of / the ground state configuration 
is given by the second line of Eq. (|36D and the DoS follows the odd asymptotics. 

This qualitative argument predicts that for asymptotically small energies, the DoS scales 
as in Eq. (pQ]). However, for larger values of the energy, rj ~ |1 — /|, a crossover to the 
characteristics of the odd / = DoS profile takes place. 

VII. SUMMARY 

In this paper, transport and spectral properties of weakly disordered quantum sublattice 
wires have been explored from a fieldtheoretical perspective. We re-derived results obtained 
previously within the DMPK transfer matrix formalism, observed a surprisingly strong sen- 
sitivity of system properties to the realization of the lead/device coupling, and found that 
conductance and DoS, at least of short systems, exhibit drastic dependence on the parity 
of the total site number in the sublattice chain. It is likely that both this phenomenon and 
the dependence of system properties on the parity of the channel number root in the same 
origin, i.e. the existence of zero energy states for effectively block off-diagonal Hamiltonians 
with rectangular (non-quadratic) block structure. Although we are not aware of an intuitive 
explanation for the channel number parity effects, this belief is supported by the observation 
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that all staggering phenomena are controlled by the same topological term 5'top in the ac- 
tion of the cr-model. From these findings one might expect that the delocalized band-centre 
behaviour exhibited by the two-dimensional sublattice model, is driven by the 2(i-analogue 
of this operator. Curiously, this is not so. In the two-dimensional field theory, the so-called 
Gade term S'cade, a two- gradient operator contribution with small and non- universal cou- 
pling constant, drives the system towards de-localization. The operator Stop does have a 
generalization to two dimensions P|, but its role has not been investigated so far. Sum- 
marizing we find that comparable phenomenology (metallic behaviour and diverging DoS) 
in one and two dimensions is described by different operators in the cr-model action. This 
indicates that some 'deeper' physical principle, not understood at present, lies beyond the 
visible structure of the field theory. 
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APPENDIX A: FIELD INTEGRAL FORMULATION 

cr-model representations of zero- and two-dimensional systems with Alll-symmetry have 
been constructed in different contexts before . That the following two sections dis- 



cuss the construction of the field theory in some detail is motivated by non-generic features 
particular to the Id-system, most notably the coupling operators and the existence of topo- 
logical structures. In the present section, we will derive a representation of the model and 
its correlations functions in terms of a supersymmetric field integral. The projection of this, 
a priori exact representation onto its low energy sector will be the subject of the subsequent 
Appendix 0. 

Consider the two correlation functions 
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C^pa'p' = (G'a/3(0+)G'«'/3'(0~)) conductaiice, 

relevant for the computation of DoS and conductance, respectively. To compute these ob- 
jects, we follow the now standard supersymmetry scheme for disordered electronic systems 
|3l| and represent the Green functions as 



where 



Z\.J] = ' 



J=0 



Z[Jl 



Z[J] = I P(^,^)e^'^(^-^+^-Sc^^^^"-^)'^, (Al) 

and we have assumed that sgnlm^; > 0. (In the opposite case, the sign of both the total 

action and the coupling operator change.) Here = {{Sa,Xa)^} and = {{Sa,Xa)} are 

two-component superfields where Sa is the complex conjugate of Sa, while Xa and Xa are 

independent Grassmann variables. The source field 

> 

J= ( 

J' 



where J^'^ = {Jat'} are ordinary matrices in site and orbital space. 

From this representation, the one-point correlation function obtains as C^]^, 



a' a 



J {Z[J]) for z = E^. However, at first sight it looks like two Gaussian field integrals 
l|) are needed to compute the two-particle correlator C^'^\ one for the Green function 
G{0^) the other for G{0~). Fortunately, this is not so, a direct consequence of the chiral 
symmetry of the Hamiltonian: the relation [as,H]+ = implies that 

G{z) = -a,Gi-z)a, 

and thus GaisiO'^) = (— )'"*"'^"'"^G*q/3(0~), where / and k are the site indices carried by the 
composite variables a and (3, respectively. In other words, the retarded and the advanced 
Green function are not independent and we can obtain the correlation function C*^^-* as 

6^ 



r^(2) _ (_\l'+k' 

'■' IBa"^ 13' a' 
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{Z[J]) 
3=0 



from the comparatively simple generating functional of the one-point function (evaluated 
for z = 0+). 



APPENDIX B: DERIVATION OF THE FIELD THEORY 



In this Appendix we derive the effective Lagrangian (|T0[) from the basic representation 



(|Al|). To keep the discussion simple, we will suppress the source- field dependence of the 
partition function in much of what follows. The final results for the correlation functions 
C-^'^-*, obtained by straightforward expansion of the action to first and second order in J are 



displayed in the final Eq. (pl| ). 

The derivation of the effective action essentially follows the standard construction 
route of field theories of disordered electronic systems. We begin by averaging the partition 
function over the Gaussian distribution of the random hopping matrices R: 



where Hq = {tii>6aa'} is the clean part of the Hamiltonian. Next, the quartic contribution 
is decoupled by means of two auxiliary fields: 

Z[0] = J I)(g,P)e-^S^ str(Q?,^,+P5^j J J^^^^^y^{z-Ho+^^J:a^^^n^ X 

where = Q ± iP, and Q and P are two-component supermatrix fields (reflecting the 
two-component matrix structure of the dyadic products ipil)) living on the non-directed links 
of our system. Both, internal structure and symmetry properties of these fields will be 
discussed momentarily. At this stage, we merely anticipate that the field configurations 
relevant to the long range behaviour of the system will be smooth. The structure of the 
action then suggests to define Qf = |(<5n_,_i + Qu-i) &s a new field variable, which now 
again lives on the sites of our system. 
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To account for the staggered even/odd structure of the theory we next 'double the unit 
cell' and define a two component field 

Here we have introduced a new counting index j = 0, . . . , L/2 enumerating the doubled unit 
cells of our system. To avoid confusion, we will systematically designate the index 0, . . . , L 
of the 'primitive' sites by /, /',... and the new index by j,j', . . .. Expressed in terms of ^ 
the functional integral assumes the form 

Z[0] = J p(g,P)e-^^S. str(Q|+P/) J ^^^^^^ ^ 

X exp (i^ (z + 2X{Q + tPas) + iir W^W^^ + a+H^'^ + a^H^A ^"j . 

Here, a± = (7i±a2 where ai,i = 1, 2, 3 are Pauli matrices acting in the space defined through 
the two-component structure of In lack of better terminology, we will refer to this space 
as the 'chiral space'. The second line of the expression above, purely off diagonal in the 
chiral space, contains the clean Hamiltonian. The explicit lattice structure of the blocks 
= {Hq'^Y is given by 

H^l> = 6uil + a)+6ii>-i{l-a), 
Hl\, = 6u,{l + a)+Sii,+i{l-a), 

where a is the staggering parameter introduced in (Q). Finally, notice that the chiral matrix 
structure of the coupling operator is given by 



3 ] 



At this stage, the superfield \& can be integrated out and we arrive at 

X exp l^-A^str \n(^z + ivr ^ W^W^^ + 2A(Q + iPa^) + a+H]^ + a^Hl^^ 
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To make further progress, we subject the functional integral to a saddle point analysis. 



are matrix-diagonal and spatially uniform. Further, we temporarily set a = W = Re 2; = 
and neglect boundary effects due to the finite extent of the system. Making the ansatz 
Q = iq ■ t., P = p ■ 1., where q and p are complex numbers and 1 is the two-dimensional 
unit matrix in superspace, a variation of the action w.r.t. Q and P then generates the set 
of equations 



where Gq = {iS + ct+Hq'^ + o--Hq^)^^\^^q and the trace extends over the two chiral compo- 
nents of the operators (but not over superspace). To evaluate the trace, we Fourier transform 
to momentum space. Defining the transform through 



(where the factor of two in the exponent serves as a mnemonic indicating that we have 
doubled the unit cell of our system) one finds that the Green function Gq is diagonal in 
momentum space with G^ik) = (i5 + (1 + e^*'^)(7+ + (1 + e~^*^)cr_^ . It is now straightfor- 
ward to verify that our equations are solved by p = and q determined through the self 
consistency equation 



particular representative of a whole manifold of solutions. To explore the morphology of 
that manifold, we notice that our restricted {z = W = 0) action is invariant under the 
transformation 



Following the standard scheme [^, we seek for saddle point configurations Q and P that 




/Ox 1/^ 
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^1 f^^^i, ^2 ^ Ti^2, 

Q + ^ Ti(g + tP)T2, Q~iP^ T^\Q - iP)T{\ 

where Ti,T2 G GL(1|1). This symmetry, the super-generahzation of the fermion symmetry 
(0), states that the model has GL(1|1) x GL(1|1) as a global invariance group. Applying 
the transformation to the diagonal stationary phase solution discussed above we find that 

—a —i^ 

where T = T1T2. Arguing in reverse, we conclude that any matrix T defines a solution of the 
mean field equations, i.e. GL(1|1) is the Goldstone mode manifold of the model. This is the 



celebrated mechanism of chiral symmetry breaking (see Ref. for review): field theory 
implementations of models with a discrete chiral symmetry on the microscopic level, possess 
continuous factor groups G x G as symmetry manifolds. (In our case G = GL(1|1).) This 
symmetry is spontaneously broken by the saddle point configurations of the model. What 
remains is a Goldstone mode isomorphic to a single factor G. 
Combining these results, we parameterize our field manifold as 

Q-iP) 2 \^ T-^U ^ 
where both T,U & GL(1|1). Here, the matrices T span the Goldstone model manifold 

whereas the U's, incompatible with the global symmetry of the model, represent massive 

modes. The next logical step in the construction of the field theory is to substitute these 

configurations back into the action and to expand in (i) energy arguments z and matrix 

elements W, (ii) long-ranged spatial fluctuations of the Goldstone modes, and (iii) massive 

modes. 



1. Goldstone mode fluctuations 



We begin with the second element of the program formulated above, the expansion of 
the action in long ranged spatial fluctuations of the Goldstone mode. Temporarily setting 
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z = W = 0, U = '1, we re-organize the 'str In' of the action according to 



fiX^T Hi'' \ ( Hl^' 

X = Nstr In = A^str In 

\H^' iX'T-^j \THl^T-^ 

where T is a slowly fluctuating fleld of Goldstone modes. Writing THq^T~^ = Hq^ + 

T[Hl^ ,T~^] and noticing that, due to the slow fluctuation of T, the commutator is small, 

we expand as 

X = Nsti {G^'T[Hl\ T-^]) - ystr {G^'T[HI\ T-^]G''T[H^\ T-']) + . . . , 

where we have deflned Go = (Go + iX)^^ and the ellipses stand for infrared irrelevant 
higher order commutator terms. One next Fourier transforms these expressions, substitutes 
the explicit momentum representation Go{k) and uses that the characteristic momentum q 
carried by the transforms T{q) is small. The subsequent integral over the 'fast momentum' 
k is most economically done by noticing that full integration over k amounts to integrating 
the characteristic phases exp(i2A;) once over the complex unit circle. This integral has a 
simple pole inside the integration contour, whose residues depend on the 'small' momentum 
q. Expanding the residues to lowest non-vanishing order in q and transforming back to 
coordinate space one obtains X = Stop + Sq, where the two contributions are displayed in 
Eq. (|lTD and a continuum limit J2j ~^ | lo dr has been taken. 

Notice that the above discussion implicitly assumed that the number of sites of our 
system is even: the operator X had a structure where each T (living at an odd site) came 
with a partner at the neighbouring even site. For a system with an odd number of 
sites L, however, there remains one uncompensated degree of freedom T((L — l)/2) at 
the terminating site. Neglecting gradients, the action due to this extra contribution reads 
'S'top,2 = A^str InT(L). In principle, this contribution is as relevant as the 'bulk' contribution 
to the topological action: it is local in space but, on the other hand, does not contain 
derivatives. Indeed, the structure of 5'top,2 is closely related to that of Stop as can be seen 
by representing the latter as a boundary action (cf. the third line of (O).) Yet in the 
majority of cases, the extra contribution Stop, 2 does not play much of a role, wherefore we 
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have largely ignored it in main analysis. E.g., for a system coupled to the outside world, 
'S'top,2 = A^str ln(T(L)) = N2ni, evaluates to a phase of no physical effect (cf. the related 
discussion around Eq. (p!3|).) There is, however, one exception to that rule, viz. the physics 
of isolated short systems with an odd number of sites discussed in section |V1 A| , where the 
presence of the extra contribution is of key relevance. 



2. Finite energies and coupling to the leads 

To obtain the action associated to finite z and WW'^, we organize the action as, 
str In 



' z + iTiWW^ + iX^T HP 







H^"" z + inWW^ + iX'^T-^ 

str In 



z + i7cWW^)T-^ \ f iX^ HP' 

+ 



{z + mWW^)T I \HP iX^ 



str [{z + t7T{WW^)2j+i)Tr'G]] + {z + t7T{WW%j)TiGf^ 



= -izTTUo Str (T + T-i) + ^ E str [iWW^)2j+iTr' + {WW^)2jTi 

where we have used that (from the saddle point equation) u = i^G]^ = i^G]] = ^ and u 
is the DoS per site evaluated for energies far away from the middle of the band. We next 
assume that in the coupling region to the leads, i.e. the region where the envelope function 
/ defined through is non-vanishing, fluctuations of the Goldstone modes are negligible. 
Application of the orthogonality relation (|) then directly leads to 

i=0,L/2 

for the contribution of the coupling term. Combining terms and taking the continuum limit, 
we finally obtain the two expressions and St of Eq. (pH]) for the contribution of energy 
and coupling operators to the effective action, respectively. 
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3. Integration over massive modes 



We finally turn to the discussion of the role played by the massive modes U. First, notice 
that due to the presence of a weight term ~ exp(— y str (Q^ + P'^)) in the action and 

AT\2 

iVstr + p^) = TVstr ((Q + iP){Q - iP)) = — str (f/^) 

8 

fluctuations of these fields are strongly inhibited. Starting out from an ansatz U = expiW, 
where W is some generator, we may thus perturbatively expand the action around W = 0. 
The actual realization of this program is cumbersome (cf. |^ for a concrete example.) Since 
the result, an extra Goldstone mode operator weakly coupled to the action, will not play 
much of a role in the present analysis we restrict ourselves to a schematic outline of the 
calculation. 

Perturbative expansion of the action in powers of W obtains an expression like 



Z[0] = y PTe""^!^! ^g-5(i)[T,vy]+5(2)[r,iy]+... 



W ^ 



where {. . .)w = J VWe-—'^'^^'^ and S^''^[T,W] denotes the expansion of the action to 
n-th order in W. The ellipses stand for contributions of higher order in W which can safely 
be neglected (due to the large overall factor ^ 1). It is straightforward to verify that 
the explicit evaluation of the operators S''-"^ [W, T] obtains contributions of the structure 
ciA^ str (14^$) and C2N str (W^W^), where $ = TdT~^ and the coupling constants ci^2 
are proportional to powers of A. Substituting these expressions back into the action and 
performing the Gaussian integration over W, we arrive at 

(str(5(^)[T, W]f)w (X Ncl str {^^) = Ncj sir (dTdT-^), 
{str {S^^\T,W]))w oc C2 [str($)]^ = C2 [str (T^T"^) " 

What can be said about structure and relevance of these operators? First, the two ex- 
pressions str ($^) and str^($) are the only invariant Goldstone mode operators with two 
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derivatives. The first of tliese is already contained in tlie action (cf. Eq. (|TT])) witli a 
coupling constant parametrically larger than the constant Ncl obtained above. Thus, the 
first of the two contributions coming from the massive mode integration is irrelevant. In 
contrast, the second contribution must be taken seriously and, after integration over spatial 
coordinates, gives the term ^cade of Eq. (PH). 

This completes our derivation of the effective action of the model. Finally, to compute 
the correlation functions C^^'^-* defined in the text, we have to add the source field J to the 
partition function, expand to first or second order and differentiate. The structure of the 
resulting expressions depends on the index configuration of the correlation functions. For 
the correlators relevant to the computation of conductance and DoS, respectively, we obtain 



Ti , / even 

I (T^-\,^ , /odd 



2 



Ti 1 I even, /' even 

^4,i2^(T'-i) ) , /even, /'odd 



(Bl) 



where the angular brackets stand for the functional average. 



APPENDIX C: GEOMETRY OF GL(1|1) 

The canonical metric on the supergroup GL(1|1) derives from the differential two-form 
UJ = — sti (dTdT^^). UJ is not a positive two-form, but its restriction to the sub-manifold 
M. C GL(1|1) defined in the text is; Ai is the maximally Riemannian subset of GL(1|1). 
Parameterizing the group manifold in terms of some coordinates, T = T[xi, . . . , X4) (half of 
which are anti-commuting), the metric two-form assumes the form 

UJ = ^Qijdxidxj, 
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where {gij} defines the metric tensor (represented in the basis {xi}.) As with non-super Rie- 
mannian manifolds, the metric tensor determines the geometry of the manifold. Specifically, 
the invariant group integral is defined through 

j dT = J Y[dxi J{xi, . . . ,0:4), 

i 

where the Jacobian J oc sdet^f. Similarly, the Laplacian has the standard structure 

where g^^ is the inverse of g, gijg^^ = 5f . 

We next wish to represent the metric tensor, the invariant measure, and the Laplacian in 
the polar decomposition, T = kak~^ defined in Eq. (|l^). Due to the manifest invariance of 
uj under transformation T — > k^TkQ^, k^ a fixed rotation matrix, it is sufficient to evaluate 
the metric tensor at A; = U. Substitution of the decomposition T = kak~^ into the two form 
then yields the covariant structure 



UJ 



str 



(d{kak-^)d{ka'^k-^\ 



str {{[dk, a] + da) {[dk, a ^]+ da 



— str 



{\^dk, a] [dk, a '^]+ dada . 



Using that 



dk 



dr] ' 



e^dx 



da 



du J \ ie^^dy ^ 

it is straightforward to verify by elementary matrix manipulation that the metric bilinear 

form assumes the form 

/I \ 



1 



V 



+4sinh2(^^ 







where g = {gij}, a vectorial structure (ix = {dx, dy, drj, dv) is understood and empty blocks 
are filled with zeros. The associated superdeterminant, 

1 



g 



[a sinh^ ( 
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x—iy 
2 



2 ■ 



From g derives the unit-normalized group integral 

/POO rliT p 

dTf{T) = dx dy dr]diyj{x,y)f{x,y,r],iy), 
J-oo Jo J 

with Jacobian 

J{x,y) = sinh-2 (~~^~^) • 
Further, the radial part of the Laplacian reads as 

i=x,y 

APPENDIX D: CONDUCTANCE OF SHORT SYSTEMS 

Consider a short sublattice system of length L < ^/{M'-y). This is the quantum dot 
regime, where the conductance is not Ohmic but rather determined by the couphng of the 
system to the leads. For a normal, non-sublattice system, g ~ M7, reflecting that each of 
the M channels contributes to transport with an efficiency set by the coupling. We here wish 
to explore to which extent this behaviour generalizes to the sublattice system. Specifically, 
three different cases will be discussed: (a) smoothened coupling where a set of sites at each 
end is connected - the type of coupling considered in the text, (b) coupling through a single 
site on either end where the two terminal sites are of the same parity, odd/odd, say. (c) 
One of the two terminating sites is even, the other odd. 

That the system is short means that the functional integral is controlled by spatially 
uniform configurations T = const.; the stiffness introduced by the gradient term is too 
strong to allow for significant fiuctuations. (More precisely, fiuctuating field configurations 
lead to relative corrections of 0{LM'y/C,) which we are not going to consider.) 

Let us begin by considering case (a). Evaluation of the functional expectation value (||) 
for a spatially uniform field configuration leads to the expression 
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(CI) 



(C2) 




We next need to do the group integral. Inspection of the exponent shows that the integration 
is dominated by configurations close to the group origin, T = U. This suggests to represent 
our T's as T = exp W and to integrate over the generators in a Gaussian approximation: 

where we have used that close to the group origin the integration measure is fiat. Doing the 
integral over the components of W one finds 

Ja) = Mil 

^ 2 

in agreement with the behaviour of a non-sublattice system. The situation in case (b) is 
not much different. Noticing that fields T sit at the odd sites of our system, while T~^'s are 
attached to even sites, we find that the conductance is expressed as 

(The disappearance of the combination T — T^^ in the pre-exponent refiects the fact that 
only single sites of different parity are coupled to the continuum. Again the integral is 
dominated by the group origin and similar reasoning as above leads to 



2 

In case (c), however, the situation is different. Owing to the identical parity of the termi- 
nating sites, the integral representation of the conductance now assumes the form 



This is an unpleasant expression: the exponential weight no longer projects onto the group 
origin, in fact does not even have a stable saddle point. To compute the conductance we 
therefore have to integrate over the full group manifold, a task that is most efficiently done 
in polar coordinates. Using Eqs. (|T^ and ( |(J1| ) and integrating out Grassmann components 
it is straightforward to verify that the radial part of the integral assumes the form 



9 



2tt (qX _ ^iy\2 



'W_oo Jo sinh2(^) 
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To understand the vanishing of the integral, notice that I dy — f y ^® transformed 
into the integral of the complex variable z = e*^ over the unit circle. One verifies that, 
regardless of the value of x, the integrand is analytic and void of singularities inside the 
integration contour. Cauchies theorem then implies vanishing of the integral. 
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